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ABSTRACT 

We use deep HST/ACS observations to calculate the star formation history (SFH) of the Cetus 
dwarf spheroidal (dSph) galaxy. Our photometry reaches below the oldest main sequence turn-offs, 
which allows us to estimate the age and duration of the main episode of star formation in Cetus. This 
is well approximated by a single episode that peaked roughly 12±0.5 Gyr ago and lasted no longer 
than about 1.9±0.5 Gyr (FWHM). Our solution also suggests that essentially no stars formed in Cetus 
during the past 8 Gyrs. This makes Cetus' SFH comparable to that of the oldest Milky Way dSphs. 
Given the current isolation of Cetus in the outer fringes of the Local Group, this implies that Cetus 
is a clear outlier in the morphology-Galactocentric distance relation that holds for the majority of 
Milky Way dwarf satellites. Our results also show that Cetus continued forming stars through z ~ 1, 
long after the Universe was reionized, and that there is no clear signature of the epoch of reionization 
in Cetus' SFH. We discuss briefly the implications of these results for dwarf galaxy evolution models. 
Finally, we present a comprehensive account of the data reduction and analysis strategy adopted for 
all galaxies targeted by the LCID (Local Cosmology from Isolated Dwarfs 15 ) project. We employ two 
different photometry codes (DAOPHOT / ALLFRAME and DOLPHOT), three different SFH recon- 
struction codes (IAC-pop/MinnlAC, MATCH, COLE), and two stellar evolution libraries (BaSTI and 
Padova/Girardi), allowing for a detailed assessment of the modeling and observational uncertainties. 
Subject headings: Local Group galaxies: individual (Cetus dSph) galaxies: evolution galaxies: pho- 
tometry Galaxy: stellar content 



1. INTRODUCTION 

A powerful method to study the mechanisms that drive 
the evolution of stellar systems is the recovery of their full 
star formation history (SFH). This can be done by cou- 
pling deep and accurate photometry, reaching the oldest 
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main sequence (MS) turn-off (TO), with the synthetic 
color-magn i tude d iagram (CMP) mo d eling technique 
dTosi et all Il991t iBertelli et al.1 Il99l iTolstov fe Sahal 



19961 lAparicio et alJ Il997t lHarris fc Zaritskvf 



Dolph in 2002; Aparicio fc GaHartll2004HCole et all 



2001] 
20071 : 



Aparicio fc Hidaigoll2009[ ). The details of the first stages 
of galaxy evolution are particularly interesting, be- 
cause they directly connect stellar populations research 
with cosmological studies. For example, the standard 
paradigm predicts that the role of re-ionization can im- 
pact the star forma tion history of small systems in a mea- 
surab le way (e.g., Ikcuchi 1986; Rces 1986; Efstathiou 
1991 iBabul fc Reesl LL991 IChTba fc Nathl H99T 
Quiim Katz. fc Efstathiou! 119961; iThoul fc Weinberg 
19961: iKepner. Babul, fc Spergell Il997t iBarkana fc Loet 



19991: Bullock et all 120001: iTassis et al.1 12003 
Ricotti fc Gnedinl 120051 : iGnedin fc Kravtsovl 120061 : 
Okamoto fc Frenkll2009ft . 

In this context, we designed a project aimed at recov- 
ering the full SFHs of six isolated galaxies in the Local 
Group (LG), namely IC 1613, Leo A, LGS 3, Phoenix, 
Cetus, and Tucana, with particular focus on the details 
of the early SFH. Isolated systems were selected because 
they are thought to have completed few, if any, orbits 
inside the LG. Therefore, compared to satellite dSphs, 
they spent most of their lifetime unperturbed, so that 
their evolution is expected not to be strongly affected by 
giant galaxies. A general summary of the goals, design 
and outcome of the LCID project can be found in Gallart 
et al. (in prep), where a comparative study of the results 
for the six galaxies is presented. 

In this paper we focus on the stellar populations of 
the Cetus dSph galaxy only. This galaxy was discov- 
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ered by I Whiting et al.l (|1999f ). by visual inspection of 
the ESO/SRC plates. From follow-up CCD observations 
they obtained a CMD to a depth of V ~ 23, clearly 
showing the upper part of the red giant branch (RGB), 
but no evidence of a young main sequence, nor of either 
red or blue supergiants typical of evolved young popu- 
lations. From the tip of the RGB (TRGB) they could 
estimate a distance modulus (m — M)o = 24.45 ± 0.15 
mag (776 ± 53 kpc) , assuming E(B — V) = 0.03 from the 
iSchleeel et aD (fl998h maps and an absolute magnitude 
M^TRGB) = -3.98 ± 0.05 mag for the tip. The color 
of the RGB stars near the tip allowed a mean metal- 
licity estimate of [Fe/H] = —1.7 dex with a spread of 
~ 0.2 dex. Deeper H ST/WFPC2 data presented by 
iSaraiedini et al.l (|2002T) confirmed the distance modulus 
of Cetus: (m-M) = 24.49 ±0.14 mag, or 790 ±50 kpc, 
assuming Mi(TRGB) = -4.05 ±0.10 mag. Based on the 
color distribution of bright RGB stars, they estimated a 
mean metallicity [Fe/H] = —1.9 dex. An interesting 
featur e of the Cetus CMD presented by ISaraiedini et al.l 
(2002) is that the horizontal branch (HB) seems to be 
more popula ted in the red than i n the blue part (see also 
the C MD bvlTolstov et all 120 00. based on VLT/FORS1 
data) . ISaraiedini et al.l (|2002f) calculated an HB index = 
—0.91 ± 0.09 16 . By using an empirical relation between 
the mean color of the HB and its morphological type, de- 
rived from globular clusters of metallicity similar to that 
of Cetus, they concluded that the Cetus HB is too red for 
its metallicity, and therefore, is to some extent affected 
by the so-called second parameter. Assuming that this 
is mostly due to age, they proposed that Cetus could be 
2-3 Gyr younger th an the old G a lactic globular clusters. 

Recent work by iLewis etldl (|2007l ) presented wide- 
field photometry from the Isaac Newton Telescope to- 
gether with Keck spectroscopic data for the brightest 
RGB stars. From CaT measurements of 70 stars, they es- 
timated a mean metallicity of [Fe/H] = — 1.9 dex, fully 
consistent with previous estimates based on photomet- 
ric data. Moreover, they determined for the first time a 
systemic heliocentric velocity of ~ 87 km s _1 , an inter- 
nal velocity dispersion of 17 km s _1 and, interestingly 
enough, little evidence of ro tation. 

The search for H I gas (Bo uchard et al.1 [2006! ) reve- 
lead three possible candidate clouds associated with Ce- 
tus, located beyond its tidal radius, but within 20 Kpc 
from its centre. Consi dering the velocity of these clouds, 
iBouchard et al.l (|2006t ) suggested that the velocity of Ce- 
tus should be of the order of -280±40fcm s" 1 to have a 
high chance of association. Cle arly, this does not agree 
with the estimates later given bv lLewis et al.l(|2007D . who 
pointed out that Cetus is devoid of gas, at least to the 
actual observational limits. 

Cetus is particularly interesting for two distinctive 
characteristics. First, the high degree of isolation, it be- 
ing at least ~ 680 kpc away from both the Galaxy and 
M31. Cetus is, together with Tucana, one of the two 
isolated dwarf spheroidals in the LG. Moreover, based 
on wid e-field INT telescope data, iMcConnachie fc Irwin! 
(2006) estimated a tidal radius of 6.6 kpc, which would 

16 The HB index is defined as (B-R)/(B+V+R), iLee et al.l 
1990), where B and R are the number of HB stars bluer and redder 
than the instability strip, and V is the total number of RR Lyrae 
stars. 



make Cetus the largest dSph in the LG. Therefore, an 
in-depth study of its stellar populations, and how they 
evolved with time, could provide important clues about 
the processes governing dwarf galaxy evolution. 

The paper is organized as follows. In Sj2]we present the 
data and the reduction strategy. $3] discusses the details 
of the derived color- magnitude diagram (CMD). In 2] 
we summarize the methods adopted to recover the SFH, 
and we discuss the effect of the photometry/library/code 
on the derived SFH. The Cetus SFH is presented in fJSJ 
where we discuss the details of the results, as well as 
the tests we performed to estimate the duration of the 
main peak of star formation and the radial gradients. 
A discussion of these results in the general context of 
galaxy evolution are presented in £j6j and our conclusions 
are summarized in fjTJ 

2. OBSERVATIONS AND DATA REDUCTION 
STRATEGY 

The data presented in this paper were c ollected with 
the ACS/WFC camera (|Ford et al.l Il998l ) aboard the 
HST. The observations were designed to obtain a signal- 
to-noise ratio f=a 10 at the magnitude level of the oldest 
M STO, Mfsi4w_ ~ ±3. Following the precepts outlined 
in IStetsonl (|1993l ) , the choice of the filters was based on 
the analysis of synthetic CMDs in the ACS bands. The 
(FA75W — F814W) pair, due to the large color baseline 
together with the large filter widths, turned out to be 
optimal for separating age differences of the order of 1-2 
Gyr in stellar populations older than 10 Gyr. 

The observations were split into two-orbit visits. One 
FA75W and one F81AW image were collected during 
each orbit, with exposure times slightly different between 
the two orbits of the same visit: 1,280 s and 1,135 s, re- 
spectively, in the first orbit, 1,300 s and 1,137 s during 
the second one. The last visit included a third orbit, with 
exposure times of 1,300 s and 1,137 s. The 27 orbits al- 
located to this galaxy were executed between August 28 
and 30, 2006. The two orbits corresponding to the ninth 
visit (images j9fz27glq, j9fz27gmq, j9fz27goq, j9fz27gqq) 
suffered from loss of the guide star. Therefore, two im- 
ages were found to be shifted by «1.5', and the other 
two were aborted after a few seconds. These data will 
not be included in the following analysis. Summarizing, 
the total exposure time devoted to Cetus and used in this 
work was 32,280 s in FA75W and 28,381 s in F8UW. 

For our observing strategy, we split each orbit into 
one F475W and one F814W exposure in order to have 
the best p ossible time sampling for short-period variable 
stars (see Ber nard et al.l 120091 ) An off-center pointing 
was chosen to sample a wide range of galactocentric radii. 
The edge of the ACS camera was shifted to the south 
by ~ 2' from the center of Cetus, in order to avoid a 
bright field star. The coordinates of the ACS point- 
ings are [RA = 00 ,t 26 m 09 s .70, DEC = -11°04'34.7"]. 
Fig. [1] shows the location of both the ACS and the 
parallel WFPC2 field [RA = 00 ,l 25 m 51 s .41, DEC = 
-11°00'32.4"]. The pointing of the latter was optimized 
to explore the outskirts of the galaxy, a nd the result- 
ing C MD has already been presented in iBernard et al.l 
(2009). Due to the much shallower photometry in the 
WFPC2 field, this will not be included in the analysis 
of the SFH. Fig. [5] shows a color drizzled, stacked im- 
age. Note the strong gradient in the number of stars 
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Fig. 1 .— 9'x 9' DSS image showing the ACS and the WFPC2 
pointings. The ellipse represents the Cetus core radius. 

when moving outward toward the south. There are also 
a few bright field stars, and a sizable number of extended 
background sources. 

In the following, we outline the homogeneous reduc- 
tion strategy adopted for all the galaxies of the LCID 
project. In particular, we stress that we performed two 
parallel and independent photom etric reduction s, based 
on the DAOPHOT/A LLFRAME (lStetsonl ll994h and the 
DOLPHOT packages polphinll2000bl) . in order to inves- 
tigate the impact of using different photometry codes to 
obtain the CMD used to derive the SFH (see j jO|) . 

2.1. DAOPHOT reduction 

After testing different approaches to drizzle the set of 
images, we opted to work on the original _FLT data. 
A custom dithering pattern was adopted for the observa- 
tions. The shifts applied were small enough (< 10 pixels) 
that the camera distortions would not cause problems 
when subsequently merging the catalogs from different 
images. The two chips of the camera have been treated 
independently, and we applied the following steps to all 
the individual images: 

• Basic reduction: we adopted the basic reduc- 
tion (bias, flat field) from the On-The-Fly-Reprocessing 
pipeline (CALACS v.4.6.1, March 13 2006). Since we 
opted for the _FLT images, we applied the pixel area 
mask to account for the different area of the sky pro- 
jected onto each pixel, depending on its position; 

• Cosmic ray rejection: we used the data quality files 
from the pipeline to mask the cosmic rays (CR). Accord- 
ing to the flag definition in the ACS Data handbook, 
pixels with the value 4096 are affected by cosmic rays. 
These pixels were flagged with a value higher than the 
saturation threshold. This allowed DAOPHOT to recog- 
nize and properly ignore them; 

• Aperture photometry: we performed the source detec- 
tion at a 3-cr level and aperture photometry, with aper- 
ture radius=3 pixels; 

• PSF modeling: careful selection of the PSF stars was 
done through an iterative procedure. Automated rou- 
tines were used to identify bright, isolated stars, with 
good shape parameters (sharpness) and small photomet- 
ric errors. PSF stars were selected in order to sample 
the whole area of the chip, to properly take into account 



the spatial variations of the PSF. The last step was a vi- 
sual inspection of all the PSF stars, for all the images, to 
reject those with peculiarities like bright neighbors, con- 
tamination by CRs, or damaged columns. On average, 
we ended up with more than 150 PSF stars per image. 

We performed many tests in order to understand which 
PSF model, among those available in DAOPHOT, would 
work best with our data. We used the same sets of stars 
for each fixed analytical PSF model, and tested different 
degrees of variation across the field, from purely analyt- 
ical and constant (variability = —1,0) to varying cu- 
bically (va = 3). We evaluated the quality of the final 
CMD, and the best results were obtained with va > 0, 
but the tightness of the CMD sequences did not vary 
appreciably when using a quadratically or cubically vari- 
able PSF. We also checked the photometric errors and 
sharpness estimators, but could not find any evident im- 
provement when increasing the degree of variation of the 
PSF. We therefore adopted a simpler and faster linear 
PSF, allowing the software to choose the analytical func- 
tion on the basis of the \ 2 estimated by the PSF routine. 
In all the images, a Moffat function with index (3 = 1.5 
was selected. 

The next step was ALLSTAR profile-fitting photom- 
etry for individual images. The resulting catalogs were 
used only to calculate the geometric transformations with 
respect to a common coordinate system, using DAO- 
MATCH and DAOMASTER. Geometrical transforma- 
tions are needed to generate a stacked image with MON- 
TAGE2. The resulting median CR-cleaned image was 
used to create the input list for ALLFRAME. Rejecting 
all the objects with \sharp\ > 0.3 was an efficient way 
to remove most of the extended objects (i.e., background 
galaxies), as well as remaining cosmic rays. After the first 
ALLFRAME was run, we used the updated and more 
precise positions and magnitudes to refine the coordinate 
transformations and to recalculate the PSFs, using the 
previously selected list of stars. A second ALLFRAME 
run was finally ran with updated geometric transforma- 
tions and PSFs. 

The calibration to the VEGAM AG system was don e 
following the prescriptions given bv lSirianni et all ([2005) , 
using the updated zero points from iMack et all <|2007l ) 
(because Cetus was observed after the temperature 
change of the camera in July 2006). Particular atten- 
tion was paid to calculate the aperture correction to the 
default aperture radius of 0.5". Aperture photometry 
in various apertures was performed on the PSF stars of 
every image, once the neighboring stars had been sub- 
tracted. The growth curve and total m agnitude of th e 
stars were obtained using DAOGROW (|Stetsonlll990f) . 
An aperture correction was calculated for every image, 
with the result that the mean correction was 0.029 ± 
0.009 mag. Individual catalogs were calibrated, and the 
final list of stars was obtained keeping all the objects 
present in at least (half+1) images in both bands. The 
final CMD is presented and discussed in £j3j 

2.2. Crowding tests with DAOPHOT 

To estimate the completeness and errors affecting our 
photometry, we adopted a standard method of artificial- 
star tests. We added ~ 350,000 stars per chip, adopting 
a regular grid where the stars were located at the vertices 
of equilateral triangles. The distance between two artifi- 
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Fig. 2. — Stacked, drizzled color image of the Cetus field. North is up and East is left. The field of view is ~3.4x3.4 square arcmin. The 
image shows a clear gradient in the number of stars when moving toward the south. Note the considerable number of background galaxies. 



cial stars was fixed in order to pack the largest number 
of stars but at the same time to ensure that the synthetic 
stars were not affecting each other. This value was cho- 
sen to be R = (radpsF + radfu + 1) = f3 pixels. In 
this way, overlapping of the wings of the artificial stars 
was minimized to avoid influencing the fit in the core of 
the PSF. With this criterion, taking into account the di- 
mension of the CCD, we needed seven iterations with ~ 
49,770 stars each. In each iteration, the grid was moved 
by a few pixels, to better sample the crowding character- 
istics of the image. 

To verify the hypothesis that the synthetic stars were 
not affecting each other, we also performed completeness 
tests injecting stars separated by (2 x radpsF + 1) = 21 
pixels. We did not find any significant difference in the 
completeness level of synthetic stars, meaning that the 
original separation was adequate. With this criterion, 
the incompleteness of the synthetic stars only depends 
on the real stars, while the additional crowding due to 
the partial overlap of the injected stars is negligible. 



The list of injected stars was selected f ollowin g the 
prescriptions introduced by iGallart et al.l ([1996). To 
fully sample the observed magnitude and color ranges, 
we adopted a synthetic CMD created with IAC-star 
(|Aparicio fc Ga llart 2004), adopting a constant star for- 
mation rate (SFR) between and 15 Gyr, with the 
metallicity uniformly distributed between Z=0.0001 and 
Z=0.005 at any age. The stars were added on the individ- 
ual images, taking into account the coordinate and mag- 
nitude shifts. The photometry was then repeated with 
the same prescriptions as for the original images. The re- 
sults are summarized in Fig. [3J where the Amag=m&gi n - 
mag otl t of the artificial stars and selected completeness 
levels are shown for the two bands, as a function of the 
input magnitude. 

2.3. DOLPHOT reduction 

DOLPHOT is a general-purpose photometry pack- 
age adapted from the WFPC2-specific package HSTphot 
(|Dolphinll2000bl ). 17 For this reduction, we used the ACS 
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Fig. 3. — Difference Amag=magi rl -mag ou t for synthetic stars. 
The two vertical lines mark the 90% and 50% completeness level. 
A total of ~ 700, 000 stars were simulated. 



module of DOLPHOT and followed the recommended 
photometry recipe provided in the manual for version 
1.0.3. The procedure is fairly automated, with the soft- 
ware locating stars, adjusting the default Tiny Tim PSFs 
based on the image, iterating until the photometry con- 
verges, and making aperture corrections. Similar to the 
procedures with DAOPHOT, we used the original _FLT 
images, corrected for the pixel-area mask and cosmic 
rays. We used a drizzled image from the HST archive 
as the reference for the coordinate transformations. The 
search for stellar sources was performed with a thresh- 
old <=2.5(T. The photometry was done fitting a model 
PSF calculated with TinyTim. DOLPHOT also auto- 
matically estimated the aperture correction to 0.5". The 
same calibration as for the DAOPHOT photometry was 
used. 

2.4. Crowding tests with DOLPHOT 

Artificial-star tests were preformed by injecting ~ 
140,000 synthetic stars in the original images. The in- 
put colors and magnitudes of the synthetic stars cover 
the complete range of the observed colors and magni- 
tudes. The stars uniformly sample the range —1 < 
M Fi75W - M$ liW < 5 mag, -7 < M Fi75W < 8 mag. 
The photom etry was performed fol lowing the procedure 
suggested bv lHoltzman et al.l (|2006f ) . The synthetic stars 
were injected one at a time into the images to avoid af- 
fecting the stellar crowding in the images. Stars were 



17 The code is publ icly 

http:/ /purcell. as. arizona. edu/dolphot/ . 
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Fig. 4.— Same as Fig. |3] but for the DOLPHOT photometry. 
A total of ~140,000 stars were simulated. 



uniformly distributed to cover the whole area of the cam- 
era. Fig. 2] shows Amag=magi„-mag otlt of the artificial 
stars together with the 50% and 90% completeness levels 
for the two bands. 

2.5. Comparison of the two photometry sets 

A direct comparison between the DAOPHOT and 
DOLPHOT photometry is presented in Fig. [5] We de- 
tect a small zero point offset for the brightest stars in 
both bands ((m Fi75W ,DAO - m Fi7bWtDO L) ~ 0.023 mag, 
for m Fi75W < 25, and (m F siAW,DAO ~ m F8 uw,DOL) ~ 
—0.010 mag, for m F suw < 24), and a trend as a func- 
tion of magnitude: moving toward fainter magnitudes, 
DAOPHOT tends to give slightly fainter results (up to 
Am_F475w — 0.05 mag, at m F 4 7 5w = 29.1, Ampsuw — 
0.07 at m F suw — 28.0). Note that since the trend with 
magnitude is similar in the two bands, the color differ- 
ence shows a residual zero point offset (~ 0.04 mag), 
but no dependence on magnitude. We performed an ex- 
haustive search for the possible cause of such small dif- 
ferences, but its origin remained elusive, and it is likely 
a result of the many differences between the approaches 
of the two photometric packages (PSF, sky determina- 
tion, aperture correction calculation). Note that similar 
trends were detected and di scussed by iHoltzman et al.1 
<|2006ft and iHill et al.l (|1998f ) when comparing different 
photometric codes applied to HST data. This reinforces 
the idea that small differences in the resulting photome- 
try naturally arise when using different reduction codes, 
at least when dealing with HST data. 

However, we highlight that, for this project, the 
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Fig. 5. — Magnitude and color difference between the DAOPHOT 
and DOLPHOT photometry, as a function of the DAOPHOT mag- 
nitude. We detect a systematic residual zero point between the two 
codes, and a small trend as a function of magnitude. At the faintest 
end, DAOPHOT gives fainter measurements than DOLPHOT. 
However, since the trend is similar in the two filters, the color 
difference showed only a small residual zero point of the order of 
0.04 mag. The plume of objects at mi?475W ~ 25.5 is made of RR 
Lyrae variable stars. 

MinnlAC/IAC-pop method presented in 14. H is not sensi- 
tive to zero-point systematics, and therefore the impact 
of small differences in the photometry on the derived SFH 
is minimized. 

3. THE COLOR-MAGNITUDE DIAGRAM 

Fig. |5] presents the final CMD from both the 
DAOPHOT and DOLPHOT photometry, calibrated to 
the VEGAMAG system. This is the deepest CMD ob- 
tained to date for this galaxy, spanning more than 8 mag, 
from the TRGB down to «1.5 mag below the oldest 
MSTO. The limiting magnitude (F814W ~ 28.8 mag) 
is similar in both diagrams. No selections were applied 
to the DAOPHOT CMD other than those already ap- 
plied to the star list, while the DOLPHOT photome- 
try was cleaned according to the sharpness (\sharp\ < 
0.1), crowding (crowd < 0.7), and photometric er- 
rors {o- mFi75W ,a mFsliw < 0.2). The DAOPHOT and 
DOLPHOT CMDs contain - 50,600 and 44,900 stars, 
respectively 

A simple qualitative analysis reveals many interesting 
features. The TO emerges, for the first time and very 
well characterized in our photometry, at mpsuw ~ 27 
mag. The morphologies of the TO and the sub-giant 
branch region suggest that Cetus is a predominantly old 
system. There is no evidence of bright main sequence 
stars indicative of recent star formation. A plume of rel- 
atively faint objects appears at 26 < ttifsmw < 28, 
< (m,Fi75W — wiF8i4w) < 0.5, and could indicate either 



a low level of star formation until ^2-3 Gyr ago, or the 
presence of a popula tion of blue stra ggler stars (BSs) (cf., 
Mapelli et aTT l2007t iMomanv et~all l2007t iMapelli et atl 



2009). We noticed that the over-density right above the 



TO, at 0.8 < (mF475W—mF8i4w) < 1, 26.5 < itifsuw < 
27.0, is mostly present in the central regions. We verified 
that, with the DAOPHOT photometry, such overdensity 
appears as well in the CMD of the synthetic stars used in 
the completenss tests and located in the most crowded 
areas of the images. Since we simulated the same CMD 
in all the regions of the camera, there is no reason to 
expect differences in the output synthetic CMDs as a 
function of the position. This indicates that the higher 
crowding in the central region is responsible for this fea- 
ture. Nevertheless, it is possible that some of the stars 
located in this region are physical binaries. 

Concerning the evolved phases of stellar evolution, the 
RGB between 20 < mpsuw < 26.5 mag appears as a 
prominent feature. The width of the RGB is similar in 
both diagrams, and suggests some spread in age and/or 
metallicity. The sudden increase in the RGB width at 
n^Fsaw < 23.2 may be explained by the superposition of 
the RGB and the asymptotic giant branch (observational 
effects like saturation of bright stars or loss of linearity 
of the camera can be safely excl uded in this regime). 

T he HB, as already no ted by ISaraiedini et al.l (2002) 
and iTolstoy et al.l (|2000f ) is dominated by stars redder 
than the RR Lyrae instability strip. However, even if 
sparsely populated, a blue HB emerges clearly from our 
photometry. It is difficult to estimate how extended it is 
on the blue side, because it merges with the sequence 
of blue candidate MS stars/BSs. We also note that 
there are a small number of objects on both the blue 
(rriF8i4W < 23 mag, (mF475W — mF8i4w) ~ 1) and red 
side of the RGB. These are most likely foreground field 
stars. 

We can gain some insight into the Cetus stellar popula- 
tions through a comparison with theoretical isochrones, 
which allows us to estimate some rough limits to the 
ages and metallicities of the stellar populations. Fig. 
[7| presents the superp osition of selected i sochr ones from 
the BaSTI database (jPietrinferni et alj [20041 ) . for the 
indicated ages and metallicities. For simplicity, we 
adopted the DAOPHOT photometry, but the conclu- 
sions are independent of the photometry set. The dis- 
tance, (to — M)q = 24.49, was est imated using the me an 
magnitude of the RR Lyrae stars ([Bernard et al.lr2~009j ) 18 . 
An attempt to use the TRGB showed that the number 
of stars near the tip is too small to derive a precise es- 
timate of the distance modulus. However, the excellent 
agreement between the TRGB of the isochrones with the 
brightest stars of the RGB in Fig. [7] shows that this dis- 
tance is very reli able. The exti n ction E(B — V) = 0.03 
was taken from ISchlegel et al.1 (|1998f L an d was trans- 
forme d to Af475\v and Afsuw following iBedin et al.l 
([2005h . 

The left panel of Fig. [7] shows the comparison with 
young isochrones. The blue edge of the blue plume can 
be well matched with a metal-poor (Z = 0.0001), rela- 

18 The value adopted in the SFH derivation, 24 .49, is slightly 
different from the final value of 24.46 ± 0.12 given in ( Bernard et al. 
2009). However, such a small difference has negligible impact on 
the derived SFH. 
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Fig. 6. — Color-magnitude diagram, calibrated in the VEGAMAG photometric system, for both the DAOPHOT and DOLPHOT 
photometry. On the right side of both panels we report the mean magnitude and color fitting errors provided by the two codes. 



tively young (2 Gyr) isochrone (red line). A similarly 
good agreement with isochrones of higher metallicity (Z 
= 0.001, 0.002, green and blue lines) is obtained only by 
reducing the age considerably (0.5 Gyr), but the lack of 
observed bright TO stars makes this interpretation un- 
likely. An increase in the age for the Z = 0.001 isochrone 
(pink line) cannot account for the observationally well 
measured blue portion of this feature. This might sug- 
gest that this is a sequence of BSs rather than a truly 
young MS (see also N5.2.ip . 

The right panel of the same figure shows the compari- 
son with isochrones of age ranging from 7 to 14 Gyr, and 
metallicities Z = 0.0001 and Z = 0.001. It is interest- 
ing to note that both the TO and RGB phases are well 
confined in this metallicity range. In particular, the blue 
edge of the RGB seems to correspond to the lower metal- 
licity, while the red edge is well constrained by a metal- 
licity 10 times higher, with negligible dependence on age. 
The situation around the TO is more complicated. At 
fixed metallicity, the oldest isochrone matches the faint 
red edge of the TO reasonably well, while the youngest 



nicely delimits the brighter and bluer one. Therefore, 
both the faint red and the bright blue edges can be prop- 
erly represented using either an older, more metal-poor, 
or a slightly younger, more metal-rich isochrone. 

However, from this simple analysis, we can not set 
strong constraints. A possible interpretation is that Ce- 
tus consists of a dominant population of age 10 Gyr, 
with a small age range and with metallicity in the range 
0.0001 < Z < 0.001. Alternatively, and equally consistent 
with this simple analysis, is a scenario with both age and 
metallicity spreads. Therefore, this simplistic analysis is 
only helpful to put loose constraints on the expected pop- 
ulations, in terms of ages and metallicities, but does not 
allow us to draw any quantitative conclusions on the age 
and metallicity distribution of the observed stars. For 
that, a more sophisticated analysis is required such as 
the one presented in g] below. 

3.1. The RGB bump 

The RGB bump is a feature common ly observed in old 
stellar systems (e.g.. iRTello et al1l2003l ). whose physical 



Monelli et al. 



i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i 

DAOPHOT (m-M) = 24.49 
.Z = 0.0001, 2 Gyr 
7 = 0.001, 0,5 Gyr 
Z = 0.002, 0.5 Gyr 
7 = 0.001, 2 Gyr • 



i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 
E(B-V) = 0,03 

.Z = G.O001, 14 Gyr 

Z = 0.001 , 10. Gyr 

Z = 0.0001, 10 Gyr 

Z = 0.001, 7 Gyr ■" 



20 



22 



5 



£ 24 



26 



28 





-0,5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3. 

m F475W ~~ m F814W m F475W ~~ m F814W 



Fig. 7. — Color-magnitude diagram, calibrated in the VEGAMAG photometric system, with overimposed some selected isochroncs from 
the BaSTI database, shifted according to the labeled distance modulus and reddening values. Age and metallicity of the isochrones are 
also labeled in the figure. 



origin is relatively well understood (|Thomaslll967t llbenl 
1968). When a low-mass star ascends the RGB, the H- 
burning shell moves outward and reaches the chemical 
discontinuity left by the convective envelope at the time 
of its greatest depth. This forces a rearrangement of the 
stellar structure and a temporary decrease in luminosity, 
before the star resumes evolving toward higher luminosi- 
ties and lower effective temperatures. From the obser- 
vational point of view, since the star crosses the same 
luminosity interval three times, it produces an accumu- 
lation of objects at one point in the luminosity function 
of the RGB. Fig. [8] shows the observed Cetus RGB lumi- 
nosity function for both sets of photometry The sample 
stars were selected in a box closely enclosing the RGB, 
in order to minimize the contribution of AGB and HB 
stars. The RGB bump clearly shows up, and fitting a 
Gaussian profile indicates a peak magnitude of tofsi4 ~ 
23.84 ±0.08, with excellent agreement between the two 
photometry sets. 



4. DERIVATION OF THE STAR FORMATION 
HISTORY OF CETUS 

In this section we present the various approaches used 
to derive the SFH of Cetus. In particular, we adopted 
the two sets of photometry previ ously described, two stel- 
lar evolution libraries, BaSTI 19 (iPietrinferni et a,l.ll2004f) 
and P adova/Girardi 20 (|Girardi et al.ll2000HMarigo et alj 
l2008h. and three differe nt SFH co des: IAC-po p 
(jAparicio fc Hidalgo! l200?)h. M ATCH (|Dolnhinl 12003) . 
and COLE (Skillman^TaD [2001 . We did not 
explore all the possible combinations of photome- 
try/library/method, but we derived different solutions 
useful for comparison purposes, and to search for possi- 
ble systematics. In particular, IAC-pop was applied to 
both sets of photometry and libraries, while MATCH and 
COLE were tested with the DOLPHOT photometry in 
combination with the Girardi library. 



http: //www. oa-teramo. inaf. it/ BASTI 



http://pleiadi.oapd. inaf.it/ 
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• a set of bolometric corrections to transform the the- 
oretical stellar evolution tracks into the ACS cam- 
era photo metric system . We a pplied the same set, 
taken from lBedin et al.l (|2005| ). to both libraries. 

• the SFR as a function of time, tp(t). We used a 
constant SFR in the age range < t < 15 Gyr; 




Fig. 8. — Luminosity function of the Cetus RGB, for both the 
DAOPHOT (solid line) and DOLPHOT (dashed) photometry. The 
vertical line shows the position of the RGB bump. 

4.1. The I AC method 

This method has been developed in recent years by 
members of the Instituto de Astrofisica de Canarias 
(A. A., S.H., C.G.), and now includes different indepen- 
dent modules. Each software pac kage has been pre- 
sente d in different papers (IAC-star: [Ap aricio fe Gallartl 
2004 MinnlAC: Hidalg o et al., in prep., IAC-pop: 
Aparicio fe Hid algo 2009), which are the reference publi- 
cations for the interested reader and include appropriate 
references 21 . In the following we give only a practical 
overview of the parameters adopted to derive the SFH of 
Cetus. 

The entire procedure can be divided into five steps. 
Note that in the following we will use "synthetic CMD" 
to designate the CMD generated with IAC-star, and we 
will call it a "model CMD" if the observational errors 
have been simulated in the synthetic CMD. 

1 ) The synthetic CMD : it is created using IAC-star 
(|Aparicio fe Gallartl l2004| ) . We calculated a synthetic 
CMD with 8,000,000 stars from the inputs described in 
the following. We have verified that, with this large num- 
ber of stars, no systematic effects in the d erived SFH 
of mo ck galaxies — such as those discussed in iNoel et al.l 
(2009) — were observed, and so no corrections such as 
those described in that paper are necessary. The re- 
quested inputs are: 

• a set of theoretical s tellar evolution models . We 
adopted the Ba STI (iPietrinferni et "all I2004D and 
Padova/Girardi (jGirardi et al.l l2000) libraries; 

21 IAC-star and IAC-pop are freely 

availa ble at http://iac-star.iac.es and 

|http: //www. iac.es / galeria/ aaj / iac-pop_eng.htm , respectively. 



• No a priori age-metallicity relation Z(t) is adopted: 
the stars of any age have metallicities uniformly 
distributed in the range 0.000K Z <0.005. We 
stress that the age and metallicity intervals are de- 
liberately selected to be wider than the ranges ex- 
pected in the solution. This is done to ensure that 
no information is lost, and to provide the code with 
enough degrees of freedom to find the best possible 
solution; 

• the in itial mass function (IMF), taken from lKroupal 
(2002). This is expressed with the formula 
N(m) dm — mT a dm, where a = 1.3 for stars 
with mass smaller than 0.5 M Q , and a — 2.3 for 
stars of higher mass. Different values of both ex- 
ponents have been tested with IC 1613 and LGS 3. 
The results, presented in Skillman et al. (in prep.), 
show that the best solutions are obtained with val- 
ues compatible with the Kroupa IMF; 

• the binary fraction, /3, and the relative mass distri- 
bution of binary stars, q 22 . We tested six values 
of the binary fraction, from 0% to 100%, in steps 
of 20%, with fixed q > 0.5; in the following analy- 
sis we use 40% and discuss the impact of different 
assumptions in the Appendix I; 

• the assumed mass loss durin g the RGB pha se fol- 
lows the empirical relation bv lReimersI ([1975), with 
efficiency r\ = 0.35. 

2) The error simulation: The code to simulate the ob- 
servational errors in the synth etic CMD, cal l ed obs ersin, 
has been developed following iGallart et al.1 (|1996h , and 
is described in Hidalgo et al., in prep. It takes into ac- 
count the incompleteness and photometric errors due to 
crowding using an empirical approach, taking the infor- 
mation from the completeness test described in £ 12 . 21 and 
^2.41 This is a fundamental step because the distribu- 
tion of stars in the observed CMD is strongly modified 
from the actual distribution due to the observational er- 
rors, particularly at the faint magnitude level near the 
old TO, which is where most of the information on the 
most ancient star formation is encoded. Fig. |9] shows 
an example of synthetic CMD as generated by IAC-star 
(upper panel), and after the dispersion according to the 
observational errors of Cetus (lower panel). 

3) Parameterization: MinnlAC is a suite of routines 
developed specifically to accomplish two main purposes: 
first, an efficient sampling of the parameter space to ob- 
tain a set of solutions with different paramctrizations of 
the CMD, and second, estimation of the average SFH 

22 This limit is set for practical purposes, since binaries with 
high er mass ratios do not a ffect the distribution of stars in the CMD 
(see IHurlev Sz Tout 1998). A 40% fraction of bina ries with this 
appro ach implies a larger total binary fraction. See Gall art et al.l 
(1999) for a discussion of this and related issues. 
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Fig. 9. — Hess diagram of the 8 X 10 6 stars synthetic CMD gener- 
ated with IAC-star (upper) and after simulating the observational 
errors (lower). 

from the different solutions that best represent the ob- 
servations. In our case, "parameterization" of the CMD 
means two things: i) defining the age and metallicity 
bins that fix the simple populations and ii) defining a 
grid of boxes covering critical evolutionary phases in the 
CMDs, where the stars of the observed and model CMD 
are counted. 

The basic age and metallicity bins defining the simple 
populations used in this work are: 
age : [1.5 2 3 ... 13 14 (15)]*10 9 years 
metallicity : [0.1 0.3 0.5 0.7 1.0 1.5 2.0]*10" 3 

Note that model stars younger than 1.5 Gyrwere not 
included because there is no evidence of their presence 
in Cetus. These numbers define the boundaries of the 
bins, not their centers. Therefore, they fix 15 x 6 = 90 
simple populations. The choice of a 1 Gyr bin size at 
all ages was adopted after testing a range of bin widths. 
It was found that smaller bins do not increase the age 
resolution; rather they increase the noise in the solution. 
The selected size of 1 Gyr is the optimal compromise to 
fully exploit the data, within the limits imposed by the 
observational errors and the intrinsic age resolution in 
the CMD at old ages. 

The parameterizat ion of the CMD relies on the con- 
cept of bundles (see lAparicio"fc Hidalgo 120091 and Fig. 
ITDj) . that is, macro- regions on the CMD that can be sub- 
divided into boxes using different appropriate samplings. 
The number of boxes in the bundle determines the weight 
that the region has for the derived SFH. This is an effi- 
cient and flexible approach for two reasons. First, it al- 




Fig. 10. — Cetus observed CMD, transformed to the absolute 
plane assuming (m-M) = 24.49 and E(B-V)=0.03, with the lo- 
cation of the four bundles superimposed. 



lows a finer sampling in those regions of the CMD where 
the models are less affected by the uncertainties in the 
input physics (MS, SGB). Second, the dimensions of the 
boxes can be increased, and the impact on the solution 
decreased, in those regions where the small number of ob- 
served stars could introduce noise due to small number 
statistics or where we are less confident in the stellar evo- 
lution predictions. We optimized the values for this work 
using four bundles, with the following box sizes. The 
smallest boxes (0.01 mag x 0.2 mag, Nl ~ 900 boxes) 
sample the lower (old) MS, and the SGB (bundle 1). 
The region of the candidate BSs (bundle 2) is mapped 
with slightly bigger boxes (0.03 mag x 0.5 mag, N2 ~ 
90 boxes), to minimize the number of boxes with zero 
or few stars. Bundle 3 samples the region between the 
BS sequence and the RGB. No bundles include the stars 
in evolved evolutionary phases, and in particular we ex- 
cluded both the RGB and the HB from our analysis. The 
reasons are many. First, the physics governing these evo- 
lutionary phases is more uncertain than that describing 
the MS, and differences between stel lar libraries are mor e 
important in these evolved phases (Gal lart et al.l f2005). 
Moreover, the details of the HB morphology also depend 
on highly unknown factors like the mass loss during the 
RGB phase. However, despite not including the RGB 
stars, we routinely adopt a bundle which is redder than 
the RGB (bundle 4 in Fig. IT0|) . In order to aid in setting 
a mild constraint on the upper limit of the metallicity 
of the RGB Cetus stars, this was drawn at the red edge 
of the RGB. It includes metal rich stars present in the 
model, but very few objects of the observed CMD. Both 
bundles 3 and 4 contain boxes of 1.5 mag x 1.0 mag (N3 
= 7, N4 — 4). Extensive tests have shown that including 
a bundle on the RGB did not improve the solution, nor 
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produce a solution significantly different. Rather, the % 2 
of the solution increases substantially and likely spurious 
populations appear in the 3D histogram, such as old and 
metal-rich stars. 

For a given set of input parameters, the simple popu- 
lations plus the grid of bundles+boxes, MinnlAC counts 
the stars in the boxes for both the observed and all the 
simple populations in the model CMD. A simple table 
with these star counts is the input information to run 
IAC-pop. However, our numerous tests disclosed that 
this approach can be significantly improved by adopting 
a series of slightly different sets of input parameters, used 
to derive many solutions and then calculating the mean 
SFH. This "dithering approach" , which is the distinctive 
feature of the code, significantly reduces the fluctuations 
associated with the adopted sampling (both in terms of 
simple populations and boxes) of each individual solution 
(Hidalgo et al., in prep.). 

In the case of the Cetus SFH, the age and metallic- 
ity bins are shifted three times, each time by an amount 
equal to 30% of the bin size, with four different configura- 
tions: i) moving the age bin toward increasing age (with 
fixed metallicity) ; ii) moving the metallicity bin toward 
increasing metallicity (with fixed age); Hi) moving both 
bins toward increasing values; and iv) moving toward 
decreasing age and increasing metallicity. These 12 dif- 
ferent sets of simple populations are used twice, shifting 
the boxes a fraction of their size across the CMD. 

Moreover, to take into account the uncertainties associ- 
ated with the distance and reddening estimates, together 
with all the hidden systematics possibly affecting the 
zero points of the photometry, MinnlAC also repeats the 
whole procedure after shifting the observed CMD (not 
the model) a number of steps in both color and magni- 
tude. The bundles are correspondingly shifted. For this 
project, we adopted an initial grid of 25 positions, shift- 
ing the observed CMD in magnitude by [—0.15, —0.075, 
0, +0.075, +0.15] mag, and in color by [-0.06, -0.03, 0, 
+0.03, +0.06] mag. Since we calculated 24 solutions for 
each node of this grid, we derived a total of 600 solutions. 

We stress that the position of the best xt m the mag- 
nitude and color grid is not intended for estimates of 
distance or reddening, since photometry zero points or 
model systematics also play a role. However, we think it 
is reasonable to take the best agreement between model 
and observations in that particular point of the {5 C oi, 
Smag) grid as pointer of the best "absolute" SFH solu- 
tion. 

4) Solving: IAC-pop (jAparicio fc Hidalgo! I2009D is a 
code designed to solve the SFH of a resolved stellar sys- 
tem. It uses only the information of the star counts in the 
different boxes across the CMD. Therefore, it calculates 
which combination of partial models gives the best rep- 
resentation of the stars in the observed CMD, using th e 
modified xl merit-function introduced bv lMigheH (1999). 
Note that IAC-pop solves the SFH using age and metal- 
licity as independent variables, without any assumption 
on the age-metallicity relation. IAC-pop was run inde- 
pendently on each different initial configuration created 
with MinnlAC. 

5) Averaging: A posteriori, MinnlAC was also used 
to read the individual solutions, and calculate the mean 
best SFH. 



4.2. The MATCH method 

The basis of the MATCH method of measuring SFHs 
(Dolphin 2002) is that it aims at minimizing the differ- 
ence between the observed CMD and synthetically gen- 
erated CMDs. The SFH that produces the best-fit syn- 
thetic CMD is the most likely SFH of the observed CMD. 
To determine the best-fit CMD, the MATCH method 
uses a Poisson maximum likelihood statistic to compare 
the observed and synthetic CMDs. For specifi c details 
on how this is implemented in MATCH, see iDolphinl 
(I2000aft . 

To appropriately compare the synthetic and observed 
CMDs, we binned them into Hess diagrams by 0.10 
mags in mp^bw and 0.05 mags in uip^w — ^ifsuw- 
For consistency, we chose parameters to create synthetic 
CMDs that closely matched those used in the IAC-pop 
code presented in this paper: 

• a power law IMF with a = 2.3 from 0.1 to 100 M ; 

• a binary fraction (3 = 0.40; 

• a distance modulus of 24.49, Af 475ty = 0.11; 

• the stellar evolution libraries of iMarigo et al.l (|2008f) . 
which are the basic Girardi models with updated AGB 
evolutionary sequences. 

The synthetic CMDs were populated with stars over a 
time range of logt = 6.6 to logi = 10.20, with a uniform 
bin size of 0.1. Further, the program was allowed to solve 
for the best fit metallicity per time bin, drawing from a 
range of [M/H] = —2.3 to 0.1, where M canonically repre- 
sents all metals in general. The depth of the photometry 
used for the SFH was equal to 50% completeness in both 
m ,F475W and mp 8 i4\Y- 

To quantify the accuracy of the resultant SFH, we 
tested for both statistical errors and modeling uncertain- 
ties. To simulate possible zero point discrepancies be- 
tween isochrones and data, we constructed SFHs adopt- 
ing small offsets in distance and extinction from the best 
fit values. The rms scatter between those solutions is 
a reasonable proxy for uncertainty in the stellar evolu- 
tion models and/or photometric zero-points. Statistical 
errors were accounted for by solving fifty random realiza- 
tions of the best fit SFH. Final error bars are calculated 
summing in quadrature the statistical and systematic er- 
rors. 

4.3. The COLE method 

The third method we applied to measure the SFH of 
Cetus is a simulated annealing algorithm applied to our 
DOLPHOT photometry and the Padova/ Girardi library. 
The co de is the same a s that applied to the galaxies 
Leo A (jCole et al.l I2007H and IC 1613 (|Skillman et al.l 
2003), and details of the method can be found in those 
references. The technique treats the task as a classic 
inverse problem by finding the SFH that maximizes the 
likelihood, based on Poisson statistics, that the observed 
distribution of stars in the binned CMD (Hess diagram) 
could have been drawn from the model. 

The code uses computed stellar evolution and atmo- 
sphere models that give the stellar colors and magni- 
tudes as a function of time for stars of a wide range 
of initial masses an d metallicities. The isochrone tables 
(M arigo et aLll2008f ) for each (age, metallicity) pair are 
shifted by the appropriate distance and reddening values, 
and transformed into a discrete color-magnitude-density 
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Fig. 11. — Summary of the different results applying different techniques to derive the SFH. The cumulative mass fraction as a function 
of time is shown. Top: Comparison of the SFH derived with the MinnlAC/IAC-pop method, using the BaSTI library and the DAOPHOT 
and DOLPHOT photometry. The DAOPHOT SFH is overall less than 1 Gyr older then the DOLPHOT one. Middle: The same solutions 
are compared to analogous ones derived using the Padova/Girardi library. These solutions are systematically younger than the BaSTI 
solutions, particularly at old epochs for the DAOPHOT photometry. Bottom: Comparison of two different SFH codes, applied to the 
DOLPHOT photometry in combination with the Padova/Girardi library. 



distributio n p by convolu tion with an initial mass func- 
tion (see iChabrierl I2003L , their Tab. 1). and the re- 
sults of the artificial star tests simulating observational 
uncertainties and incompleteness. The effects of binary 
stars on the observations are simulated by adopting the 
binary star frequency and mass ratio distribution from 
the solar neighborhood, namely 35% of stars are single, 
46% are parameterized as "wide" binaries, i.e., the sec- 
ondary mass is uncorrelated with the primary mass, and 
is drawn from the same IMF as the primary, and finally 
19% of stars are parameterized as "close" binaries, i.e., 
the probability distribution function of the mass ratio is 
flat. The progeny of dynamical interactions between bina- 
ries (e.g., blue stragglers) are not modelled. The history 
of the galaxy is then divided into discrete age bins with 
approximately constant logarithmic spacing in order to 
take advantage of the increasing time resolution allowed 
by the data at the bright (young) end. The star forma- 



tion history is then modelled as the linear combination of 
V>(t,Z) that best matches the data by summing the values 
of g at each point in the Hess diagram. The best fit is 
found via a simulated annealing technique that converges 
slowly on the maximum-likelihood solution without be- 
coming trapped in local minima (such as those produced 
by age-metallicity degeneracy). 

In our solution we considered only the DOLPHOT pho- 
tometry. The critical distance and reddening parameters 
were held fixed at (m-M) = 24.49 and and E(B-V) 
= 0.03. We used 10 time bins ranging from 6.60 < 
log(age/Gyr) < 10.18. Isochrones evenly spaced by 0.2 
dex from —2.3 < [M / H] < —0.9 were used, considering 
only scaled-solar abundances but with no further con- 
straints on the age-metallicity relation. The CMD was 
discretized into bins of dimension 0.10x0.20 mag, and it 
is the density distribution in this Hess diagram that was 
fit to the models. The error bars are determined by per- 
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turbing each component of the best-fit SFH in turn and 
re-solving for a new best fit with the perturbed compo- 
nent held fixed. The size of the perturbation is increased 
until the new best fit is no longer within la of the global 
best fit. 

4.4. Comparison of the different solutions 

Fig. [11] summarizes the comparison of the solutions 
calculated using different sets of photometry, stellar evo- 
lution libraries, and SFH codes. 

The upper panel shows the cumulative mass fraction, 
as a function of time, based on the BaSTI library + IAC 
method applied to the two sets of photometry. The solu- 
tion based on the DOLPHOT photometry is marginally 
younger than the DAOPHOT one. It shows slightly lower 
star formation at older epochs, and a steeper increase 
around 12 Gyr ago. In either case Cetus completely 
stopped forming stars ~ 8 Gyr ago. 

The middle panel compares the same two curves with 
analogous ones obtained using the Girardi library. The 
effect of changing the library introduces a small system- 
atic effect; both solutions with the Girardi library shift 
more star formation to slightly (< 1 Gyr) younger ages 
compared to the solutions using the BaSTI library. This 
effect seems more evident at the oldest epochs in the case 
of the DAOPHOT solutions. 

The bottom panel shows the results of the three differ- 
ent SFH reconstruction methods, applied to the same 
photometric catalog and with the same stellar evolu- 
tion library. Since for the MATCH and COLE method 
we present the result of the best individual solution, 
and not the average of many as in the case of the IAC 
method, the large dots highlight the centers of the age 
bins adopted. Note the good agreement between the 
MATCH and COLE with the IAC solution, in spite of 
the higher time resolution of the latter. 

For simplicity, we will now focus our analysis on the 
SFHs obtained with the IAC method and the BaSTI stel- 
lar evolution library. At this time, we favor using the 
BaSTI stellar evolution library because the input physics 
has been more recently updated , in comparison to the G i- 
rardi stellar evolution library (jPietrinferni et all [2004). 
Similarly, the IAC method will be used preferentially be- 
cause it was tailored to the needs of this project and 
affords a greater degree of flexibility for the analysis (Hi- 
dalgo et al., in prep.). 

5. THE SFH OF CETUS 

In this section we present the details of the Cetus SFH. 
First, in M5 . H we describe the analysis we did to explore 
the possible systematics affecting the solution. We then 
describe our final solution in 



5.1. Approaching the best solution 

5.1.1. Exploring the zero points systematics 

Fig. [12] shows the grid of color and magnitude shifts 
adopted to derive the IAC-pop solution. The plus signs 
mark the positions of the 25 initial points. In each of 
these 25 positions we calculated the mean xt> averaging 
the xt °f the 24 individual solutions. This allows us to 
study how the xt varies as a function of the magnitude 
and color shifts. In particular, the minimum averaged 
xt min identifies the position corresponding to the best 
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Fig. 12. — Summary of the solutions derived as a function of 
the shifts in color and magnitude, <5 co ; and 5 m ag, applied to the 
observed CMD. In each position of the grid we calculated 24 solu- 
tions. Crosses mark the initial grid, common to both photometry 
sets. Filled circles and open squares show the sub-grid around the 
initial minimum, for the DAOPHOT and DOLPHOT photometry 
set, respectively. The two curves show the 1-<t confidence area, 
defined using the dispersion of the 24 individual solution around 
the mean xt- 



solution, which is calculated by MinnlAC averaging the 
corresponding 24 individual solutions. 

Because the initial color and magnitude shifts are rela- 
tively large, we calculated more solutions on a finer grid 
with 14 positions around the identified minimum xt (Fig. 
[T2l full circles and open squares for the DAOPHOT and 
DOLPHOT solution, respectively). Therefore, we ob- 
tained a total of (25+14)*24=936 solutions. This compu- 
tationally expensive task was made possible by the Co n- 
dor workload management system (iThain et al.l l2005h 23 
available at the Institute dc Astrofisica de Canarias. In 
the case of the DAOPHOT photometry, the minimum 
was confirmed at the position (5 co i ; S mag ) = (0.0 ; 
0.075), xl,min = i- 95 ' while in the case of the DOLPHOT 
photometry the minimum is located at (S co i ; 5 mag ) — 
(+0.03, 0.075), xl, mi „ = 2.88. The two curves around 
the minimum position of Fig. [12] represent the 1-er con- 
fidence area around the minimum, calculated as the dis- 
persion a round the mean of the 24 x 2 of the individual 
solutions (jAparicio k, Hi dalgo 2009|) . Note that the color 
difference between the xt min 01 the two photometry sets 
is fully consistent with, and apparently compensates for, 
the color shift between the two photometry sets (~ 0.04 
mag) . This indicates that the flexibility of our method is 
able to deal with subtle photometric calibration system- 
atics, and possible stellar evolution model, distance, and 
reddening uncertainties. 

Fig. [13] shows a direct comparison between the 24 in- 
dividual DAOPHOT and DOLPHOT solutions at the 
xt min position. The horizontal lines show the size of 
the age bins, and do not represent error bars. It is note- 
worthy that the spread of the 24 solutions, derived using 
different input samplings, is quite small, suggesting the 
robustness of the mean SFH. Note also the good agree- 
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Fig. 13.— Comparison of the 24 individual DAOPHOT and 
DOLPHOT solutions at the best xt position, from both sets of pho- 
tometry. Note the stability of the solution against small changes 
of the input bins adopted. Moreover, the plot puts into evidence 
the similarity of the solutions obtained from the DAOPHOT and 
DOLPHOT photometry. Note that the horizontal lines arc not the 
error bars, but represent the size of the age bin. 
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Fig. 14. — Results of a test designed to investigate the signifi- 
cance of the oldest t/>(t) in our solutions. Dashed line: input SFH, 
adopted as our Cetus solution, but removing all stars older than 
13.5 Gyr. Solid line: recovered SFH when using a model CMD 
that contains stars as old as 15 Gyr. A low level of star formation 
is recovered for ages older than 13.5 Gyr, when there should be 



ment between the solutions from the two photometry 
sets. The general trend is the same, with the main peak 
of star formation occurring within ~ 0.5 Gyr and, most 
importantly, characterized by the same duration. 

5.1.2. The age of the oldest stars 

Another important aspect to be addressed is the oc- 
currence of star formation at epochs older than 14 Gyr. 
As discussed in §4.11 we did not impose any strong con- 
straints on the ages of the oldest stars in the model used 
to derive the SFH of Cetus. That is, we let the code 
search for the best solution using model stars as old 
as 15 Gyr, significantly older than the age of the Uni- 
verse commonly accepted after the WMAP experiment 
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Fig. 15. — Comparison between the Cetus ip(t) obtained using a 
model CMD with constant SFR up to 15 (dashed black line) and 
and up 13.5 Gyr (solid grey line). 




([Bennett et alJl2003f >. Interestingly enough, our solution 
shows very low star formation rate for epochs older than 
14 Gyr, with a steady increase toward a peak at ~ 12 
Gyr. 25 

To further investigate the significance of deriving ages 
of stars older than 14 Gyr, we performed the following 
test. Using the derived SFH, we calculated a synthetic 
CMD, from which we removed all the stars older than 
13.5 Gyr. We simulated the observational errors and 
then derived the SFH of this new CMD identically as for 
Cetus. Fig. [141 presents the input and recovered ij)(t). 
The latter shows some level of star formation at epochs 
older than 13.5 Gyr, when there should be none. This 
means that the combination of observational errors and 
the uncertainties associated with the SFH algorithms is 
responsible for the smearing of the main peak of star 
formation at the oldest epochs. 

Therefore, we decided to adopt the external constraint 
on the age of the Universe coming from the cosmic mi- 
crowave background experiments, and we set the upper 
limit to the model stars' age to 13.5 Gyr. This choice 
is also supported by the fact that current stellar evolu- 
tion models give Milky Way globular cluster ages in very 
good agreement with this estimated age of the Universe 
(jMarm-Franch et al.ll2009h . The effect of this constraint 
is shown in Fig. [15j The main effect of limiting the age 
of the oldest stars in the synthetic reference CMD is that 
the peak of star formation is higher (since total star for- 
mation must be conserved). However, there is no effect 
on the age or the end of the main episode of star forma- 
tion, nor on the star formation at epochs younger than 
10 Gyr. With both models, we calculate that Cetus com- 
pletely stopped forming stars ~ 8 Gyr ago. This test was 
repeated with both the DAOPHOT and DOLPHOT so- 
lutions, leading to the same conclusion. Thus, we adopt 
this constraint on the SFH for the rest of our analysis. 

24 It is also worth noting that none of the stellar evolution li- 
braries implemented in IAC-star accounts for the occurrence of 
atomic diffusion. Atomic diffusion has the effect of reducing the 
evolutionary lifetimes during the central H-burning stage of low- 
mass stars, resulting in a reduction of the stellar age - at the oldes t 
ages - by about 0.7 Gyr IjCastellani et al. 1997; Cassisi et al.H999T ). 
Therefore, in the present analysis, an age estimate of ~ 14 Gyr 
would correspond to an age of about 13.3 Gyr if models account- 
ing for atomic diffusion were to be adopted. Both are in very good 
agreement with current determination of the age of the Universe. 
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Fig. 16. — 3-D representation of the Cetus SFH, derived with 
the DAOPHOT (top) and the DOLPHOT (bottom) photometry. 
Both show the average of 24 solutions calculated at the minimum 
Xv in the 8mag-8col grid. The X„ m i n ls 2.18 and 2.92 for the 
DAOPHOT and DOLPHOT solutions, respectively. 



5.2. The final SFH of Cetus 

The solutions for each set of photometry, calculated as 
the average of the 24 solutions at the Xvmin> are re P re ~ 
sented in Fig. [TC] where a 3-D smoothed histogram sum- 
marizes the main features. This is a useful represe ntation 
givin g an overall view of the SFH of a system (Hodge 
1989), including the star formation rate as a function of 
time, the metallicity distribution function and the age- 
metallicity relation. It shows graphically how the old- 
est stars are more metal-poor, and how chemical enrich- 
ment proceeds toward younger ages. The two projections 
ip(z) and ip(t) are the metallicity distribution integrated 
over time (blue line) and the age distribution integrated 
over metallicity (red line). Errors have been calculated 
for tp(t) and 4>( z ) as the dispersion of the 24 individual 
solut ions used in the average (see lAparicio" fc Hidalgol 
2009). In particular, the ip( z ), shows a peak at the 
most metal poor regime, and a sharp decrease toward 
the more metal-rich tail. Only a negligible fraction of 



TABLE 1 

Integrated quantities derived for the Cetus dSph. 





DAOPHOT 


DOLPHOT 


fip(t)dt [10 6 M Q ] 

< ip(t) > [10~ 7 M Q yr- 1 pc- 2 ] 

< age > [10 10 yr] 

< \Fe/H] > [10" 4 dexl 


(1.81±0.05) 
(1.17±0.02) 
(1.14±0.02) 
(4.23±0.54) 


(1.91±0.07) 
(1.23±0.03) 
(1.12±0.02) 
(3.64±0.53) 



stars ever formed have metallicities higher than Z=10 . 

The overall agreement between the DAOPHOT and 
the DOLPHOT solutions, illustrated in Fig. [13] and in 
Fig. \W\ suggests that we can safely adopt the average of 
the two solutions for the SFH of the Cetus dSph galaxy, 
and the difference of the two as an indication of external 

24 

errors . 

Fig. [T7] summarizes the main results by comparing the 
star formation rate as a function of time (top), the age- 
metallicity relation (middle), and the cumulative mass 
function (bottom), for both sets of photometry with the 
adopted average calculated as explained in Hidalgo et 
al., in prep. The thick dashed and continuous lines rep- 
resent the DAOPHOT and DOLPHOT solutions, while 
the thin lines are the corresponding error bars, calcu- 
lated as the l-cr dispersion of the 24 solutions at xtmin- 
This is a statistically meaningful way of defining the un- 
certainities introduced by the SFH recovery meth od, as 
extensively discussed in lAparicio fc Hidalgo! ((2009). The 
thick cyan line is the average of the solutions obtained 
for the two photometries. Fig. [TH] and Fig. [T7] provide 
two complementary ways of presenting the derived SFH 
of Cetus. From these figures we see that Cetus is mostly 
an old, metal-poor stellar system, with the vast majority 
of stars older than 8 Gyr and more metal-poor than Z = 
0.001. Tab. [^summarizes the main integrated quantities 
derived for both photometry sets. 

The vertical lines in the top panel indicate the redshift 
z = 15 and z — 6, which mark the epoch when reioniza- 
tion was under way, the latter considered to be the epoch 
when the Universe was fully reionized (jBouwens et al.1 
120071 and references therein). The vertical lines in the 
middle panel of Fig. [Tfl mark the epoch when 10, 50 and 
90% of the stellar mass was formed. They indicate that 
Cetus experienced the bulk of its star formation before 
redshift « 2, with a peak at z ~ 4. 

In the following, we focus on two important aspects of 
the Cetus SFH: the sequence of candidate BSs and the 
actual duration of the main episode of star formation. 

5.2.1. Recent star formation or blue stragglers 

In we gave circumstantial evidence, based on the 
comparison with isochrones, that the objects bluer and 
brighter than the old MSTO are BSs, rather than truly 
young stars. More evidence of this comes from the analy- 
sis of the SFH presented in [16] The plot shows an event 
of very low SFR that produced a small percentage of 

24 From Fig. 1111 middle panel, it can be concluded that the 
adopted photometry has a larger influence on the SFH than the 
stellar evolution models used. The same is true with respect to 
assumptions on binaries (see Fig. 1231 top panel), or the IMF (Skill- 
man et al. in preparation). Therefore, using the difference in the 
SFH calculated from the two photometry sets as an indication of 
external errors seems justified. 
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Fig. 17. — Summary of the Cetus SFH results. The three panels represent, from top to bottom, the ip(t), the age-metallicity relation 
and the cumulative mass fraction. We show the results with the BaSTI library for both the DAOPHOT (dashed line) and the DOLPHOT 
(continuous) photometry, and the average of the two (thick cyan line). The dashed part of the cyan line representing the average of the 
two sets of photometry in the middle panel indicates the presence of BSs rather than a low level of intermediate-age and young SF. The 
thin lines represent t he error bars, while t he vertical lines in the upper panel mark redshift z = 15 and 6, that is the epochs corresponding 
to the reionization. (Bouwcns ct al. 2007). The dotted-dashed lines in the central panel mark the epochs when 10%, 50% and 90% of the 
mass was formed. 



metal-poor but relatively young (~ 2-4 Gyr old) stars. 
This population clearly does not follow the general age- 
metallicity relation. Fig. fTS] shows that the position of 
these stars in the solution CMD is fully consistent with 
the position of the blue plume on the CMD. This fea- 
ture is systematically present in all our solutions, and 
does not depend on the sampling, library, or photometry 
adopted. Given the age-metallicity relation derived for 
Cetus, there is no reason to expect a relatively young but 
very metal-poor population - this would require some- 
thing like an episode of late infall of very metal poor 
gas. This interpretation is nicely supported by the age- 
metallicity relation presented in Fig. [T7] (middle panel), 
which shows a continuous increase of metallicity with 
time, until a peak occurs 8 Gyr ago. After that, we ob- 
serve a sharp decline followed by a renewed increase. If 
we accept that blue stragglers are found in all or nearly 



all sufficiently populous open and globular star clusters 
as well as in the Milky Way field halo population, then 
there is no conflict between our data and the assertion 
that no star formation occurred in Cetus in the last 8 
Gyr; the most natural explanation for these stars is that 
they are blue stragglers belonging to the old, metal-poor 
population. Note that the contribution of BS to the to- 
tal estimated mass of Cetus is < 3%. Therefore, the 
presence of fainter BSs contaminating the TO and MS 
regions is expected to have negligble impact on the de- 
rived SFH. An in-depth discussion of the candidate BSs 
in both dSph galaxies of the LCID sample will be pre- 
sented in a forthcoming paper (Monelli et al., in prep.) 

5.2.2. The short first episode of SFH in Cetus 

The description presented in ^corresponds to the SFH 
of Cetus as derived using the discussed algorithms. It is 
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Fig. 18. — Syntetic CMD corresponding to the SFH of Cetus. 
The young metal-poor stars are highlighted with bigger points. 
Their position on the CMD is fully consistent with the candidate 
BSs sequence. 



clear that this solution is the result of the convolution 
of the actual Cetus SFH with some smearing produced 
by the observational errors, and the limitation s of the 
SFH recovery method (|Aparicio fc Hidalgoll2009t) . In the 
following, we will try to further constrain the features 
of the actual Cetus SFH through some tests with mock 
stellar populations. 

To investigate the actual duration of the dominant 
episode of star formation, and to assess the age resolution 
at the oldest ages, we performed the following test. We 
created mock galaxies characterized by fixed metallicity 
and short episodes of star formation around the mean 
value of 12.0 Gyr (see Fig. [19|). Since a Gaussian pro- 
file fit to the Cetus ip(t) yields an estimate of a = 1.53 
Gyr, we used three ij)(t) modeled by a Gaussian profile 
with a = 0.5, 1.0, and 1.5 Gyr. We simulated the ob- 
servational errors in the mock stellar populations, and 
recovered their SFH using the same prescriptions as for 
the real data. For simplicity, only the completeness tests 
of the DAOPHOT photometry were adopted to simulate 
the observational errors in the mock CMD. The number 
of stars in the mock populations was such that the total 
number of stars in the bundles was comparable to the 
case of the real galaxy. 

The first important result of this test is that the age of 
the peak of the recovered mock population is systemati- 
cally well recovered, within 0.2 Gyr. Thus, if the initial 
episode of star formation in Cetus is well represented by 
a Gaussian profile, then we are confident in our ability to 
determine the time of the peak star formation rate. Sim- 
ilar te sts have been also discussed in lAparicio &; H idalgo 
(2009) and Hidalgo et al., in prep. The second finding is 
that the duration of the recovered ij>(t) is wider than the 
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Fig. 19. — Comparison of the ip(t) derived for Cetus with that 
retrieved from three mock stellar populations. These were com- 
puted with a Gaussian profile ip(t) with <r =0.5, 1.0, and 1.5 Gyr 
peaked at 12 Gyr, as shown in the small top-right panel. The re- 
covered ip(t) is wider than the input one. This test indicates that 
the intrinsic duration of the main star formation episode in Cetus 
can be confined between the a =0.5 and a =1 Gyr. All the curves 
presented are normalized to their total area. 



input one. In particular, we measured the HWHM and 
calculated the corresponding a, which are equal to 1.44, 
1.61, and 1.95 Gyr for the three episodes of increasing du- 
ration. The increase from the input a is a consequence 
of the blurring effect of the observational errors and the 
limited age resolution at old ages. 

The comparison of the duration of the main peak in 
the Cetus and in the mock stellar populations suggests 
that the estimated duration for Cetus is an upper limit. 
Moreover, we note that this value can be further con- 
strained by the a = 0.5 and the a — 1.0 Gyr Gaussian 
profile tests, with a result close to a — 0.8 Gyr (interpo- 
lated value). Moreover, we can expect that the limited 
age resolution broadens both the older and the younger 
sides of the ip(t), around its maximum. This is supported 
by Fig. HH where we can compare the youngest epoch of 
star formation in the three models and in their recovered 
profiles. The comparison reveals that the output ip(t) is 
more extended to younger ages than the input ones, and 
that this effect is less important for increasing Ui n . This 
suggests that at least part of the star formation between 
10 and 8 Gyr ago derived in Cetus is not real. How- 
ever, we stress that the differences between the actual 
SFH and the recovered one affect only the duration of 
the star formation solution and the height of the main 
peak, but they do not alter significantly the age of the 
peak itself. 

Finally, Fig. [30] presents a comparison between the 
observed DAOPHOT CMD and a synthetic CMD corre- 
sponding to the best solution (left and cetral panels). We 
find a satisfactory agreement in all the main evolution- 
ary phases, including the RGB phase which has not been 
used to derive the solution. As far as it concerns the dif- 
ferences in the HB morphology, it is worth noting that we 
have not tried to properly reproduce this observational 
feature. Due to the strong dependence of the HB mor- 
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phology on the mass loss efficiency along the RGB, it is 
clear that a better agreement between the observational 
data and the synthetic CMD could be obtained with a 
better knowledge of mass loss efficiency and its spread. 
The right panel shows the residual Hess diagram, which 
supports the good agreement between the observed and 
the solution CMD. 

5.3. Are we seeing the effects of reionization in the 
Cetus SFH? 

One of the main drivers of the LCID project is to find 
out whether cosmic reionization left a measurable im- 
print on the SFHs of nearby, isolated dwarf galaxies. The 
Universe is t hought to have been completely reionized be- 
fore z > 6 (|Bouwens et al .11200 71. corresponding to T ~ 
12.8 Gyr ago. The comparison of these values with the 
Cetus tp(t) in Fig. [T7] (see upper panel), shows that the 
vast majority of the measured star formation in Cetus 
occurred at epochs more recent than 12.8 Gyr ago (see 
also Tab. [2]). This is also true when the actual Cetus 
SFH is considered, as discussed in section [5.2.21 There- 
fore, we conclude that there is no evident coincidence 
between the epoch when reionization was complete and 
the time in which the star formation in Cetus ended. 

However, one might ask the question of whether ac- 
ceptable solutions could be found that don't show star 
formation more recent than 12.8 Gyr ago, i.e., that are 
compatible with the hypothesis that reionization con- 
tributed to the supression of star formation in Cetus. 
With this question in mind, we performed two additional 
tests, summarized in Fig. 1211 

• First, we calculated a solution SFH of Cetus costrain- 
ing the model CMD to stars of ages older than 12.8 Gyr. 
This is not the typical approach for the IAC method, 
but nontheless the code looks for the absolute minimum 
in the allowed parameter space, and provides the best 
possible solution (grey dashed line). However, we found 
that this solution is clearly at odds with the derived best 
SFH of Cetus, and not compatible with it. Also, note 
the significantly higher xZ> which increases to 7.80, and 
the large residuals (Fig. [211 central panel). The latter, 
calculated with respect to the observed CMD identically 
as in Fig. 1201 indicate important systematics in the solu- 
tion. The excess of stars in the solution CMD in the red 
part of the TO (lightest bins), indicates the presence of 
too many stars in this region. This is in agreement with 
the overall older population, with the main peak at ~ 13 
Gyr ago. 

• We recovered the SFH of a mock population, simi- 
larly to the tests shown in Fig. [T9l created with constant 
SFR in the age range between 13 and 13.5 Gyr, and zero 
elsewhere. The solution shown in Fig. [5TJ discloses that 
the peak is recovered at the correct age, in the oldest 
possible age bin (grey solid line). This suggests that an 
actual SFH where most of stars were born before the end 
of the reionization epoch is not compatible with the so- 
lution we derived for Cetus. Moreover, the general trend 
in the residual plot (Fig. [2TJ right panel), is similar to 
the previous case, with an excess of red TO stars. 

These two tests, together with the ones shown in Fig. 
[T9l suggest that, on the basis of the available data cou- 
pled with the SFH recovery code, we can reject the hy- 
pothesis that the majority of the star formation occurred 
in Cetus in epochs older than z ~ 6. 



TABLE 2 

Mass percentiles formed in Cetus as a 
function of redshift and look-back time. 



Mass % 


redshift 


Look-Back Time (Gyr) 


10 


8.8 


13.1 


20 


5.7 


12.7 


30 


1.5 


12.3 


40 


3.7 


12.0 


50 


3.2 


11.6 


60 


2.8 


11.3 


70 


2.4 


10.9 


80 


2.0 


10.4 


90 


1.7 


9.7 



5.4. Summary 

The various pieces of information analysed in the £j5] 
can be sumarized as follows. The star formation of the 
observed area of Cetus shows a well-defined peak that oc- 
curred ~ 12 Gyr ago. By resolving the star formation of 
mock stellar populations, we showed that the actual du- 
ration of the main event is probably overestimated, due 
to the smearing effect of the observational errors. This 
does not affect the age of the main peak, but only its du- 
ration, at both older and younger epochs. Therefore, the 
last epoch of star formation, estimated at 8 Gyr ago, is 
probably actually somewhat older. The various evidence 
indicating that the blue plume of objects identified in the 
CMD consists of a population of BS, supports the idea 
that no star formation occurred in Cetus in the last 8-9 
Gyr. In Tab. [2] we present the epoch when various frac- 
tions of the total mass were formed, once we excluded 
the contribution of the BS population. However, the ac- 
tual values of the percentiles reported could be slighlty 
different due to the smearing of the solution with respect 
of the actual, underlying SFH (see H5.2.2[) . 

The derived SFH is at odds with the conclusion by 
ISarajedini et al.l (|2002| ) , who suggested that Cetus might 
be 2-3 Gyr younger that the old Galactic globular clus- 
ters. The stability of the main peak around the age of 
12 Gyr, togethe r with up-to-date estimate s of the globu- 
lar cluster ages arm- Franch et al. 2009), suggest that 
such a difference can be ruled out. 

6. DISCUSSION 

The SFH of Cetus derived in this paper shows that 
this dSph was able to form the vast majority of its 
stars coincident with and after the epoch of reioniza- 
tion. Furthermore, there are no signs in the SFH sug- 
gestive of a link between a characteristic time in the 
SFH and the time of reionization. This may be due 
to the fact that Cet us, with a veloci ty dispersion of 
a = 17 ± 2 km/s dLewis et all |2007[) is o ne of the 
brightest (M v 10.1. IWhiting et al.|[l999h . and pre- 
sumably one of the most massive dSphs in the Local 
Group. Its large mass might therefore have shielded 
its baryonic content from the effects of reionization 
(iBarkana fc Loerj|1999t lCjjraro et"aT1l2001l: [Tassis et al l 
1 20031: iRicotti fc Gne din 20051: iGnedin fc Kravtsovl [20061 : 
Oka moto fc Frenkll2009[ ). The latter reference, in partic- 
ular, argues that all classical LG dwarfs originate in the 
few surviving satellites which, at the time of reioniza- 
tion, had potential wells deeper than the threshold value 
of ~ 12 km/s. If reionization by the UV background can- 
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Fig. 20. — Comparison of the observed (left) and solution (centre) CMDs. The latter was derived extracting random stars from the 
model CMD used to derive the solution, in such a way that each simple population contributes proportionally to the calculated SFR. The 
rigth panels shows the residual Hess diagram, shown in unit of Poisson errors. The Figure disclose an overall general agreement between 
the observed and the solution CMDs, particularly in the whole TO region used to calculate the SFH. The discrepancy at the faintest 
magnitude is due to the limit magnitude of the model CMD. Note that such discrepancy does not affect the portion of the CMD included 
in the bundles, which are repeated in the left panel. The differences in the HB morphology are expected, since small differences in the 
adopted mass-loss along the RGB have a big impact on the HB morphology. 



not be related to the early truncation of star formation 
in Cetus, then other circ umstances must be a t play . Ac- 
cording to the model by lOkamoto fc Frenkl (|2009l ). the 
sudden truncation in the SFH of dwarfs like Cetus arises 
from having their gas reservoirs stripped during first in- 
fall into the potential wells of one of the massive LG spi- 
rals (M31 or MW). Another possibility frequently raised 
in the literature is that the end of star formation reflects 
the blow-out of gas driven by supernova explosions. We 

examine each of these two possibilities below. 

The radial velocity of Cetus measured by iLewis et al.l 
(2007) implies that this galaxy is receding with respect 
to the Local Group barycenter and, under reasonable as- 
sumptions about its proper motion and the mass distri- 
bution in the Local Group, that it is close to the apocen- 
ter of its orbit. Therefore, it is possible that Cetus is in a 
highly radial orbit and that it has gone through pericen- 
ter at least once. If this was the case, then it is possible 
that the truncation of the star format ion occurred due to 
the effect of tides and r am-pressure (jMaver et al.ll200U: 
lOkamoto fc Fr cnk 2009) during peric entric passage. The 



mild evidence of rotation found by ILewis et al.l (|2007l ) 



may be interpreted to imply that the progenitor of Ce- 
tus was a disk-like dwarf that underwent an incomplete 
transformation into a spheroidal syste m as a result of the 
strong tides operating at pericenter (jMaver et alj|200fl : 



lLokas et al1l2009D . The models by lLokas et al.l (|2009f) 
predict different final configurations after one passage 
for a disky dwarf. Depending on the geometry of the 
encounter, either a triaxial, rotating system or a prolate 
object with little rot ation may result. Alt ernatively, it 
has been proposed by D'O nghia et all (|2009l ) that a grav- 
itational resonance mechanism, coupling the rotation of 
the accreting disky dwarf with its orbital motion, is able 
to turn it into a dSph in a single encounter. However, 
this requires a nearly prograde encounter. The latter 
condition may not be common during the hierarchical as- 
sembly of the galaxy, but may still have occurred in just 
a few cases (prograde accretions are seen in cosmological 
hydrodynamical simulations) which would be enough to 
explain the small subset of distant dSphs. A firm mea- 
sure of the shape of Cetus might thus be an important 
step in the context of constraining its formation models. 
Independently of the mechanism driving the encounter, 
if Cetus has indeed been through the pericenter of its 
orbit, it is possible that it owes its present-day isolation 
in the LG to the fact that it gained orbital energy in the 
process and is now on a weakly bound orbit that takes 
it well beyond the normal confin es of the LG. This pos- 
sibility has been hi ghlighted by I Sales et ail ([2007D and 
iLudlow et al.l (|2009t ). who discuss the possibility that the 
odd orbital motion of Cetus might have originated in a 
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Fig. 21. — The left panel compares the i/>(t) of Cetus (black dasked line) to (i) the one of a mock population, created with stars in the 
range between 13 and 13.5 Gyr (solid grey line), and (ii) a solution derived for Cetus, but imposing the constraint to use only stars older the 
12.8 Gyr (dashed grey line). The central and right panels show the residual Hess diagram of these two solutions, calculated identically as in 
Fig. 1201 The residual CMDs show important systematics suggesting that these solutions are not compatible with the derived best solution 
of Cetus, thus supporting that no strong correlation exists between the epoch of the rcionization and the epoch of the star formation in 
this galaxy. 



multiple-body interaction during pericenter. 

On the other hand, models of the SFH in dwarfs driven 
by feedback are also able to reproduce the main prop- 
erties of dSph galaxies without invok i ng ga laxy inter- 
actions. In particular, iCarraro et al.l (|2001l) find that 
the key parameter is the dwarf's central density, while 
ISawala et al.l (120091) argue that, as in previous models 
such as lMac Low fc Ferraral f|1999T ) , low-mass systems are 
those most likely to end up as dSph galaxies. However, 
despite some differences all models of gas blow-out by su- 
pernova explosions find the mechanism to be inefficient 
above 10 7 M Q . While ultra- faint dwarfs may well be in 
the right mass regime, this is a problem for all classic 
dSphs since their initial masses were likel y at least tw o 
orders of magnitude above this threshold (|Maverl f2009). 
Note also, that models of dSph formation that are based 
only on galaxy mass and feedback are unable to repro- 
duce the morphology-density relationship. 

When compared to the dSph satellites of the Milky 
Way, the stellar content of Cetus shows striking similar- 
ities with the Draco dSph. lAparicio et al.l ()2001 ). on the 
basis of a ground-based CMD of depth similar to those 
obtained by LCID for more distant dwarfs, concluded 
that most of its star formation occurred more than 10 
Gyr ago, and that Draco stopped forming stars ~ 8 Gyr 
ago. They also infer a less intense episode ~3-2 Gyr ago, 
which is reminiscent of the population in Cetus that we 
associate with blue stragglers. Therefore, Cetus's SFH 



looks more similar to that of other nearby dSphs with 
mostly old stellar populations, such as Draco, Ursa Mi- 
nor, Leo I, Sculptor or Sextans. Their SFHs differ from 
those of other dSphs like Carina, Fornax or Leo II, where 
there is clear evidence for a substantial intermediate-age 
population. The striking observation that Cetus and 
Tucana (the other isolated LG dSph) do not follow the 
morphology-density relation of Milky Way dSph galaxies 
observed in the LG today appears to be simply an obser- 
vational coincidence. Their radial velocities indicate that 
it is possible that they were much closer to the dominant 
LG spirals when their SF ceased, indicating that their 
star formation might have been truncated by the same 
mechanism responsible for the formation of the other LG 
dSphs. 

7. SUMMARY AND CONCLUSIONS 

We used HST/ACS observations to derive the SFH of 
the isolated Cetus dSph. This effort is part of the LCID 
project, whose main goal is to investigate the SFH of iso- 
lated dwarf galaxies in the LG. To asses the robustness of 
our conclusions, the LCID methodology uses and com- 
pares the results of two photometric codes, two stellar 
evolution libraries, and three SFH-estimation algorithms. 
This allowed us to check the stability of our solution, and 
to give a robust estimate of the external errors affecting 
it (see also Hidalgo et al., in prep.). The use of differ- 
ent stellar evolution libraries and/or SFH codes does not 
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change the global picture we derived. We found that the 
BaSTI library gives solutions systematically older (less 
than 1 Gyr) than the Padova/Girardi one. The cumula- 
tive mass fraction showed in Fig. 1111 comparing the dif- 
ferent solutions shows a good general agreement within 
the error bars. Finally, the good consistency of the SFH 
solutions derived with the DAOPHOT and DOLPHOT 
photometry sets warranted to adopt the average of the 
two solutions as the final Cetus SFH. 

Our main result is that Cetus is an old and metal-poor 
system. We find evidence for Cetus stars as old as those 
in the oldest populations in the Milky Way, and com- 
parable i n age to the age of the Universe inferred from 
WMAP (|Bennett et al.ll2003t) . Star formation decreased 
after the earliest episode which likely peaked about ^12 
Gyr ago. We find no convincing evidence for a popula- 
tion of stars younger than ~ 8 Gyr old. 

A number of tests aimed at recovering the SFH of mock 
stellar populations yield a estimate of our precision for 
determining the peak of star formation of the order of ^ 
0.5 Gyr for a given stellar evolution library. Therefore, 
our estimate that the Cetus peak tp(t) occurred ^12 
Gyr ago, is a robust result. Our tests also disclosed that 
the measured duration of the star formation is artificially 
increased by the observational errors. We could estimate 
that the intrinsic extent is compatible with an episode 
approximated by a Gaussian profile with a = 0.8 Gyr (or 
1.9 Gyr FWHM), as opposed to the measured a = 1.53 
Gyr. This likely moves the end of the star formation at 
older epochs. 

The age-metallicity relation shows a mild but steady 
increase in metallicity with time, up to Z — 0.001. The 
presence of bright, blue stars in the upper main sequence 



formally implies a sprinkling of stars younger than ~ 8 
Gyr old. These stars, however, are inferred to have 
mctallicities below Z — 0.001, and thus would be clear 
outliers in the monotonic age-metallicity relation. An al- 
ternative interpretation, which we favor, is that these are 
not younger main-sequence stars but rather a population 
of blue stragglers with the same typical metallicity as the 
bulk of the population in Cetus. 

In conclusion, we do not find any direct correlation 
between the epoch of the reionization and the epoch of 
the star formation in Cetus. In particular, the reioniza- 
tion does not seem to have been an efficient mechanism 
to stop the star formation in this galaxy. Alternative 
mechanisms, such as gas stripping and tidal interactions 
due to a close encounter in the innermost regions of the 
LG, seem favored with respect of internal processes like 
supernovae feedback. 
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APPENDIX 

The effect of binary stars 

To study the impact of different assumptions on the binary content of the model CMD, we used the IAC method 
to the identically derive the SFH with six different synthetic CMDs, varying the amount of binary stars from 0% to 
100%, in steps of 20%. In all cases, we adopted a mass ratio between the two components q > 0.5, and a flat mass 
distribution for the secondary. Fig. [5U summarizes the \1 min as a function of the binary fraction for five out of the 
six galaxies of the LCID sample, based on the DAOPHOT photometry: Cetus, Tucana, LGS 3, IC 1613, and Leo A. 
The five galaxies show a common general trend: the xt decreases for increasing number of binaries, and reaches a 
plateau and possibly a minimum around 40% - 60%. We conclude that we cannot put strong constraints on the binary 
content of these galaxies on the basis of these data, although a very low binary content seems to be disfavored. 

The impact of different amount of binaries on ij;(t) is shown in Fig. [23l where we plot the cumulative mass fraction 
and the ip(t) for the six different assumptions of the binary stars applied to the DAOPHOT photometry of Cetus. 
We detect a mild trend, the solution getting slightly younger for increasing amount of binaries. This is shown also in 
the bottom panel of the same figure, where the ip(t) tends to be more extended to younger ages for increasing binary 
fraction. However, the main features are not strongly affected by the assumptions made on the characteristics of the 
binary star population. In particular, the age of the main peak of star formation is stable within ~ 0.6 Gyr. 

Moreover, we show that the position of the minimum in the xt plot as a function of S co i — 5 mag , changes as a function 
of the assumed binary content. As a general trend, we found that for increasing binary content, the minimum xt 
tends to migrate toward redder color and shorter distance, as summarized in Tab. [3] in the case of Cetus. This can 
be easily understood in terms of how the increasing number of binaries changes the distribution of stars in the model 
CMD. The higher the binaries content, the redder and more luminous the MS and MSTO. 
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